On the Entangling Power of Quantum Evolutions 
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We analyze the entangling capabilities of unitary transformations U acting on a bipartite di x di- 
dimensional quantum system. To this aim we introduce an entangling power measure e(U) given 
by the mean linear entropy produced acting with [Zona given distribution of pure product states. 
This measure admits a natural interpretation in terms of quantum operations. For a uniform 
distribution explicit analytical results are obtained using group-theoretic arguments. The behavior 
of the features of e(U) as the subsystem dimensions d 1 and di are varied is studied both analytically 
and numerically. The two-qubit case d\ = d% — 2 is argued to be peculiar. 

PACS numbers: 03.67.Lx, 03.65.Fd 



From the beginning it has been argued that entangle- 
ment is one of the crucial ingredients that allows Quan- 
tum Information processing M to outperform, for cer- 
tain tasks, any classically operating device. In this sense 
entanglement represents a uniquely quantum resource 
whose production is a sort of elementary prerequisite for 
any Quantum Computation (QC). Such a basic task is 
accomplished by unitary transformations U i.e., quantum 
evolutions acting on the state-space of the multi-partite 
system that describe non-trivial interactions between the 
degrees of freedom of the different subsystems. Even 
though almost all the unitaries satisfy this latter require- 
ment it is quite natural to ask how different IPs are 
efficient, according to some criterion to be specified, as 
cntanglers, and then by using such a criterion to analyze 
the full manifold of bi-partite quantum evolutions. 

In this paper we address this issue by introducing over 
the space of bi-partite unitaries a measure for their en- 
tangling power. This is done by considering how much 
entanglement is produced by U on the average acting on 
a given distribution of uncntanglcd quantum states. The 
kind of situation we have in mind is a procedure for en- 
tanglement production in which one randomly generates 
product states (the " cheap" resource ) according to some 
probability distribution p and then applies the transfor- 
mation U. The average entanglement obtained with the 
above scheme will be our measure e p (U) of the quantum 
evolution U. 

It is important to stress that these [/'s can represent 
different objects, both from the logical and physical point 
of view. Some prototypical instances are given by: a) A 
quantum computation using a pair of quantum registers. 
Here the entangling power measure will quantify how the 
computation U is efficient in making the first (say mem- 
ory) and the second (say computational) registers entan- 
gled. This kind of entanglement, that represents mutual 
information between the two registers, has been recently 
proved to play a role in QC viewed as a communication 



process || . b) The global evolution of a system plus its 
environment. In this case e p (U) measures the decohering 
power of the system-environment coupled evolution U. 
Engineering weak decoherence then amounts to design 
an optimal U with respect to the criterion of minimal 
entangling power, c) A single two-subsystem e.g., two- 
qubits, gate in a quantum-network. Now the entangling 
U are the two-qubit gates needed to get universal QC 0. 

To formalize our setting let us consider a bipartite 
quantum system with state space H = Hi <£> H 2 where 
dimWi = di (i = 1, 2) and U G U{H) = U{di d 2 ). If E is 
an entanglement measure over H we define the entangling 
power of U (with respect to E) as 



,(U) :=£([/ 1 Vi><8|^2» 
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where the bar denotes the average over to all the product 
states \ipi}® \ip2)- distributed according some probability 
density pi^pi,^) over the manifold of product states. 

We shall use as entanglement measure of \^>) G H the 
linear entropy 



:= l-tr lP 2 , p:=tr 2 |tf>(*|. 



(2) 



This quantity measures the purity of the reduced density 
matrix p, it can be regarded as a kind of "linearized" 
version of the von Neumann entropy S(p) = — trplnp, 
which is known to provide the essentially unique mea- 
sure of entanglement for bi-partite pure quantum states. 
One has that < E(\ty)) < 1 - l/d where the lower 
(upper) bound is reached iff is a product state (max- 
imally entangled). The measure (|^) has, with respect to 
S(p), the definite advantage of being a polynomial in \ip). 

Now we introduce some notations. We shall denote by 
Tij, = 1, ... ,4) the transposition between the i-the 
and the j-th factor of H® 2 := H 2 = (C dl ® C da ) <g> (C rfl ® 
C d2 ). Notice that T 12 and T 34 are well defined elements 
of U(TL® 2 ) only when d\ = d%, in this latter case such 
operators will be referred to as swaps. Moreover - when 
Hi = Hj - one defines the projectors := 2 _1 (l±Ty) 
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over the totally symmetric (antisymmetric) subspaces of 
Hi ® Ttj, the latter being thought of as embedded in 
TL® 2 . The space End {TL® 2 ) is endowed with the Hilbert- 
Schmidt scalar product < A, B >:= tr(A> B). Finally 
with S(TL) we shall denote the space of density matrices 
over TL. 

Proposition The entangling power ([[]) is given by 

e v (U) = 2tv[U® 2 n p U f ® 2 P u ], (3) 

where fl p := f dfx^x, ^ 2 )(|V>i)(V>i| ® |^2)(^2|)® 2 € 
S(TL® 2 ) and denotes the measure over the product 
state manifold induced by the probability distribution 

Proof. Let us observe that Eq. (j|) can be writ- 
ten in a linear form using the identity tr [(A ® B)T] = 
tt(AB) where T is the swap. Then J5(|*» = 1 - 
tr (|*)(*|® 2 T i3 ). Form this remark and the definition 
(|l|) it follows immediately (^) □ 

The result above express e p (U) as the expectation 
value over Q p of the positive operator 2 W® 2 P^ J7® 2 . 
This latter operator can be viewed as the effect asso- 
ciated to the completely positive (CP) [0 map $u on 
S(TL® 2 ) given by ^yid i-> 2 P" [/® 2 Q f/t® 2 P". This 
remark allows us to interpret the entangling power (Q) as 
probability of success of a two-party (A and B) quantum 
protocol (see Fig. 1). Suppose A (B) owns spaces TLx 
and H3 (TL2 and TL4) 

a) A and B generate pairs of states ® IV'a) accord- 
ing the distribution probability p(ipi, 1P2) [&p is prepared] 
b) Apply to each member of the pair the joint transfor- 
mation U [action of U® 2 } c) Perform a projective mea- 
surement of "\Z2Pi3- 
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FIG. 1. Scheme of the two-party protocol realizing the 
operation &u- 

Eq. (||) nicely displays several properties required 
for any entangling measure for bi-partite unitary evo- 
lutions, i) ep(C7i <g> U 2 U) = e p (U)(Ui € U{d t )). 
Indeed from the C/(di)-invariance of Pj^ one finds 
e p (Ui®U2 U)=2tv [U® 2 tt p W® 2 {U X ®U2)^ 2 P13 {Ux® 
U 2 )® 2 } = 2tT[U® 2 n p W® 2 (U 1 )\ 3 P 13 (U 1 ) 13 ] = e p (U). 



Where (Ux)i3 := Ux g> 1 ® Ux ® 1. ii) When di = d 2 , 
by denoting with T the transposition between the two 
factors of TL, one has e p (TU) = e p {U). This stems from 
T® 2 P~T® 2 = p- 4 . This leaves e p (U) unchanged, in- 
deed this label change amounts simply to the replace- 
ment tri «-> tr2 in Eq. (0). Since, for pure states, the 
two reduced density matrices are isospectral the linear 
entropy is unchanged. Moreover if £l p is swap invariant 
i.e., p(ipi,ip2) = p(?p2,i>x) onc also has e p {UT) = e p (U). 
iii) One has e p (l) = 0. This simply because Vt p P^ = 0. 
Indeed, form the definition ([|) one has f2 p Ti3 = fl p . iv) 
From the previous remarks it follows that the entangling 
power is constant along the orbits in U(TL) of the left ac- 
tion of the subgroup of the bi- local operations Ui ® E/2. In 
particular e vanishes on all the elements of such a group. 
In the symmetric case dx — d 2 the group is extended by 
the swap T. 

Different distributions p(tpx, ^2) would result in very 
different e(U). An extreme example of this obvious re- 
mark is provided by transformations U that simply per- 
mutes elements of a given basis \i) ® \j) of TL. If p is 
supported just on this basis the associated e p (U) van- 
ishes identically, while we shall show later for a different 
probability distribution that Vs can even be maximally 
entangling. Another example is given in the context of 
the case b) mentioned in the introduction. Suppose U 
admits a decoherence-free subspace C C Hx [§ the if p 
is, for any jf/^), supported in C then again e p {U) = 0. 

From now on we focus on the case in which p is the 
uniform distribution pq. With this term we refer to the 
unique U{dx) x [/((^-invariant probability distribution 
i.e., p(ipx,ip2) = p{Ux^px, U 2 Tp2)- When all the product 
state are considered to be equally easy to be prepared, 
this latter assumption on p is quite natural from the phys- 
ical point of view ||. Moreover, in view of its symmetry, 
the uniform p will result in a great computational simpli- 
fication that will allow for an explicit analytical evalua- 
tion of the average over the product-state manifold that 
appears in Eq. ([[]). 

Let us begin by proving an easy group-theoretic 
Lemma that will play an essential technical role in the 
following. 

Lemma n po = 4 C dl C d . 2 P+ P+ , CJ 1 := d (d + 1). 

Proof. Since the uniform distribution factorizes we can 
consider separately the average u>xs with respect to 
( on the first and the third factor of TL® 2 ) and the one 
LJ24 with respect to \ip2) ( on the second and the fourth 
factor of TL® 2 ) then fL, = uixz ^24- Let first observe that 
in view of definition (0) one has that f2 po is supported in 
Pjtj P24TL® 2 i.e., fl p is symmetric under the exchange of 
the first (second) and the third (fourth) factor. Moreover 
since the uniform distribution is U(dx) x U{d 2 ) invariant 
one has [U® 2 , ojxs] — 0, VZ7i S U(dx), and analogously 
for W24- Since the U® 2 's act on the totally symmetric sub- 
space irreducibly, it follows from the above commutation 
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relation and the Schur Lemma that u>x3 = 2 C P^. 
The normalization constant is found by the condition 
tro;i3 = 1. Reasoning in the same way for w 2 4 one gets 
the desired result. □ 
Proposition 1 The entangling power (Q), with respect 
the uniform distribution, is given by 



-Pa 



(U) = i-c dl c d2 Y, i*(U) 



a=0,l 



I a (U) = t(a)+ < U® 2 (T 1+Q , 3+Q ) C/t® 2 , T 13 >, (4) 

where t(a) := trTi+ Qi 3+ a . 

Proof. It is just a calculation. Insert fl po = 
C dl C d2 (1 + T i3 ) (1 + T 24 ) in Eq. (g). Notice that one 
has tr Ti 3 = dx d%, trT 24 = d\ d 2 . □ 

From the relations [Ti +Q , 3+Q , (Ux <g> £/ 2 )® 2 ] = 0, (a = 
0, 1) it follows that both the functions Iq and /i are in- 
variant under the two-sided i.e., left and right, action of 
bi- local unitaries e.g., h(U) — Ix(U Ux®U 2 ). Moreover in 
the symmetric case dx = d 2 it follows from the above that 
the entangling power (jij) can be written in a manifestly 
swap invariant form e pa (U) = 1 — C\ E;=o ^d(T l U), 
where I d {U) := d 3 + <U® 2 , T 13 U® 2 T 13 > . 
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FIG. 2. Distribution of the probability density q(e) for 
d = 2,34. 

The entangling power e defines a random variable over 
U(TL) if the latter, endowed with the Haar measure, is 
considered as a probability space. Therefore it makes 
sense to consider the associated density of probability dis- 
tribution q(e). Moreover, since the the manifold of uni- 
tary transformations over Ti. is compact, the obviously 
continuous mapping U t— ► e(U) must achieves extrema, 
in particular 3U G IA{TL): e(U) — maxy e(U). Such max- 
imally entangling U's will be referred to as optimal. In 
Fig. 2 are reported the g(e)'s obtained numerically for 
the cases d x d with d — 2,3, 4. While in the cases d > 2 
the function q(e) vanishes on both the lower and the up- 
per sides of the allowed range of e, it it remarkable that 
the two-qubit case d = 2 shows a peculiar feature: q(e) 
is a monotonic function of e. This implies that most of 



the two-qubit gates U correspond to nearly optimal ones. 
Moreover, as will be discussed later in details, the entan- 
gling power of optimal U's does not correspond to an up- 
per bound that is instead reached by all the other cases 
for d. In this sense the prototypical quantum information 
case of two qubits is quite singular. A first very natural 
question is how on the average an operator is entangling 
i.e., the mean of the q(e). 

Proposition 2 The average of the entangling power 
e Po (U) over U(d\ d 2 ) is given by 



e P0 (U) 



u (dx - 1) (da - 1) 



(5) 



d x d 2 + 1 

Proof. To prove Eq. (g) we first notice that, in view 
of definition (g), to compute the mean of the entangling 
power amounts to compute the average of the entangle- 
ment measure over all the states |$) 6 W ( not just over 
the product states). The trace of the square of the re- 
duced density matrix of |>F) is given by tr (|*)(*|® 2 T 13 ). 
Now we take the average with respect using again 
Lemma [dx d 2 {dx d 2 + l)]- 1 tr[(l + T 13 T 24 ) T 13 ] = [dxd 2 2 + 
d 2 } [d 1 d 2 {d 1 d 2 + l)]- 1 . This expression inserted in the 
definition of the entanglement measure proves Eq. (g). 
□ 

For proving bounds on the entangling power (g) it is 
useful to consider one of the states, say \tp 2 ), of the input 
product as fixed. In this case a pair of CP-maps associ- 
ated with U are naturally defined. Indeed one has (ex- 
plicit dependence on U and \ip 2 ) is omitted) $: S(TCx) > 
S(Hx):p •-» E,ti A iP A t and S(H 2 ):p ^ 

Et=i A ipA\ where Af.Hx *— ► Hx and Aj-.Hx i— > T~L 2 are 
given by A t := (j\ U \ip 2 ), (j = 1, . . . , d 2 ) 

Ai ■= EjliUX^i. (* = !,•••,*)• Therefore one 
can also define, for fixed \ip 2 ) £ TL 2 the (partial) en- 

ip 

tangling power of U as e p (<I>) := E(<& \t/j)(ip\) . No- 
tice that the equation above can also be written is 
the form (JsJ) with a special choice for p(ipt,ip2) i-e., 
with Q p = JX^XIViX^il <8 \^ 2 }(iP 2 \)® 2 . The defi- 
nition of e p (<I>) of course makes sense for general CP- 
maps, in this case the expression for e p ($) analogous to 
Eq. (g) is given by e p ($) = 2 tr [$® 2 (uj p ) P u ] where 
Lo p :=fdpW\iP)^\® 2 €S(Hf 2 ). 

Proposition 3 The entangling power of the CP-map 
Q:S(Hi) i— * S(TLx) with respect to the uniform distribu- 
tion is given by 



e P0 ($) = l-C dl (tr 2 X 2 + trxX 2 ) 



(6) 



where X := Y*L X A i A and * ■= £$=i h. A \- 

Proof. One has that trj p 2 is given by 

ES^itri^il^O^ilA^j ■ IViXV'il At), this last ex- 
pression can be rewritten as Eii=i tr i [(Aj Ai) ® 
(A\ Aj) \%pi){ipi\® 2 ]. Now we perform the average with 
respect \ipt); using Lemma u> Po — C dl (1 + T). Using 
again the identity tr [(A ®B)T]= tr (A B) one gets 
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^triP 2 ^ 1 = E |tn(4^)| 2 + tn (J^Ai A}) 2 (7) 
i,i=i 1= 

[Notice that the two terms in the equation above corre- 
sponds to the I Q 's in Eq. It is then straightforward 
algebra to check that the first term in the equation above 
can be written as tr2 X 2 □ 

We now provide bounds on the entangling power. We 
assume that d% < d%. 

Proposition 4 For any U EU (TL) one has 

S e„ (t ,> £ *^* ,8, 

Proof. The lower bound is obvious in view of the 
definition (|J), it is achieved by all the unitaries ob- 
tained composing bi-local transformations of U(d) x U(d) 
with the swap. Let us consider first the operator X 
in Eq. (||). By denoting with the CP-map associ- 
ated with the Aj's one has &(t/di) = Xjd\ and then 
1 — tri <f>(l/<ii) 2 < 1 — 1/di (general bound on linear en- 
tropy) it follows that d\ < tri X 2 . This latter inequality 
provides a bound on the second term of Eq. (0). Rea- 
soning in the same way with the operator X and the as- 
sociated CP-map one finds for the first term of (||) the 
lower bound d\jd2- Putting these two results together, 
and in view of the assumption d\ < di , inverting d\ with 
c?2 one gets the desired result (||) □ 

Another issue is to understand whether the upper 
bound in Eq. (JsJ) is achieved by an optimal unitary trans- 
formation U. As we shall show in the following the answer 
seems to be affirmative for d > 2, Let us stress that this 
is not obvious at all in that the upper bound (|^) has been 
obtained by providing separate bounds on the two terms 
appearing in e Po (U). Notice that the proof of Proposition 
3 allows us to state the condition on U in order to sat- 
urate the bound (|J) as : for any initial state \ip) 6 Hi 
the associated CP-maps $ and (depending both on 
and U) must be unital i.e., they map totally mixed 
states onto totally mixed states. It might well be that 
no £/'s yields unitality for both CP-maps at once. This 
in fact turns out numerically to be the case for d = 2, 
in which one has that the optimal i7's are such that [see 
Fig. 2] e Po (U) = 2/9 < 1/3. An optimal operators for 
qubit is given (not surprisingly) by the controlled-not 
U := |0)(0| <g> 1 + |1)(1| <g> X, where X := |0)(1| + |1)(0|. 
More interestingly the operators providing a natural d- 
dimensional generalization of the controlled-not are in 
general not optimal. This is shown by the following cal- 
culation. 

Let us consider, for di = da, U — Yla=i \ a )( a \ ®U a , 
where the |q)'s are a d-dimensional orthonormal basis 
and the U^s are unitaries which, without any loss of gen- 
erality can be taken to be orthogonal with respect to the 
Hilbert-Schmidt scalar product i.e., < U a , Up >= d5 Ut p. 
By using Eq. (0) is easy to prove that for these unitaries 



one has e PQ (U) = d(d— l)/(d+l) 2 , that is d/(d+l) times 
smaller than the bound (||) and, for d > 2 even smaller. 

We also performed numerical maximization of e Po (U) 
trying to find exact expressions for the optimal unitary 
transformations. For d± — c?2 = d = odd we have 
found that the following "classical" unitary transforma- 
tion (which only permutes the d 2 bases states) is optimal 
and reaches the bound: U \i) ® \j) = \i + j) ® \i — j) 
where i, j = ... d — 1 and the sums are mod d. [Notice 
that the above expression for even d does not define a 
permutation of the basis of H..] Similarly a more com- 
plicated construction gives an optimal permutation that 
achieves the bound for the case d\ = c?2 = d — 4n. Thus 
for equal dimensions the only case that remains to be 
solved is d\ = c?2 = d = An + 2, e.g. 6x6. We have also 
found unitary transformations satisfying the bound for a 
very asymmetric case, namely d^ — n-m with n, m > d\. 
This last example is of the type of controlled unitary op- 
eration from the larger to the smaller system. The pre- 
vious constructions for equal dimensions can be viewed 
as a concatenation of two such controlled operations, in 
which the control role is played alternatively by one of 
the subsystems. It may also be interesting that for the 
cases 2 x odd the bound can be shown not to be reached 
by permutations. 

Let us finally discuss briefly the numerical evidences. 
First, even for dimensions 2x3 the bound appears not 
to be reached, rather for optimal C/'s we get the value 
1/3 (instead of 3/8). For all other cases that we have 
checked the bound seems to be reached, namely for 2x4 
up to 2 x 7, and 3 x 4 up to 3 x 6. In conclusion one 
might conjecture that the only cases where the optimal 
transformations do not reach the bound (||) are 2x2 and 
2x3§. 

Conclusions. In this paper we introduced a measure 
for the entangling power e(U) of unitary transformations 
U acting on the state-space TL of a bi-partite d\ x c?2 quan- 
tum system. In terms of this measure we moved a first 
step towards the analysis of the manifold of bi-partite 
unitary transformations. We analyzed the induced prob- 
ability distribution q(e) over U(7i) as d\ and c?2 varies, 
and we found an analytical form of optimal transforma- 
tions for some cases. Although we believe that both the 
questions addressed and the approach we adopted are 
quite natural and physically motivated the role, if any, 
that the entangling power will play in Quantum Infor- 
mation theory is still an issue for future work. 

The authors thank M. Rasetti and J. Pachos for use- 
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IST-Q-ACTA. 
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